Density and spin response of a strongly-interacting Fermi gas 
in the attractive and quasi-repulsive regime 
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Recent experimental advances in ultra-cold Fermi gases allow for exploring response functions 
under different dynamical conditions. In particular, the issue of obtaining a "quasi-repulsive" regime 
starting from a Fermi gas with an attractive inter-particle interaction while avoiding the formation 
of the two-body bound state is currently debated. Here, we provide a calculation of the density and 
spin response for a wide range of temperature and coupling both in the attractive and quasi-repulsive 
regime, whereby the system is assumed to evolve non-adiabatically toward the "upper branch" of 
the Fermi gas. A comparison is made with the available experimental data for these two quantities. 

PACS numbers: 03.75.Ss,03.75.Hh,74.40.-n,74.20.-z 



Ultra-cold Fermi gases represent testing systems for re- 
solving many open issues in condensed and nuclear mat- 
ter. A key feature of these systems is that the inter- 
particle interaction can be varied with unprecedented 
flexibility through the use of Fano-Feshbach resonances 
from the weak- to the strong-coupling limits, which cor- 
respond to the presence of correlated and truly bound 
pairs, in the order. Recent experimental advances have 
also made it possible to achieve an accurate control of 
the temperature, in such a way that the temperature de- 
pendence of several physical quantities can be explored. 

In particular, due to the diluteness condition of an 
ultra-cold Fermi gas, the temperature interval that can 
be explored ranges from about 5% of the Fermi tempera- 
ture Tp up to several times Tp. Such a wide temperature 
range allows alternative theoretical approaches and the 
corresponding results to be tested, from the temperature 
regime T 3> Tp where the first few virial corrections to 
the free Fermi gas are relevant, down to the tempera- 
ture region T k. T c where the interplay of thermal and 
quantum fluctuations signals the presence of a supcrfluid 
phase that develops at the critical temperature T c . 

In this context, a recent experiment [l[ has reported 
values for the compressibility and spin susceptibility of a 
unitary Fermi gas over a wide temperature range (0.2 % 
T/Tp < 10), setting a benchmark for theoretical calcu- 
lations that address the (static limits of the) density and 
spin correlation functions. 

Here, we present theoretical results for the compress- 
ibility and spin susceptibility, obtained within linear- 
response theory by a diagrammatic approach built on 
the t-matrix approximation and its variations. The dia- 
grams selected for the calculation bear on familiar con- 
tributions in condensed matter, namely, the density-of- 
states (DOS), Maki-Thompson (MT), and Aslamazov- 
Larkin (AL) diagrams Q that are depicted in FigfT] 

In the temperature range T c < T < 5Tp relevant to 
compare with the experimental data, we shall obtain the 
compressibility as the static limit of the density correla- 
tion function, by adding to the DOS diagram of FigQJa) 



the MT diagram of FigQJb) > the AL diagrams of FigsOJc) 
andQJd), plus the whole series of FigJTJc) which is built 
on these AL diagrams This is because for the com- 
pressibility it is essential to take into account the residual 
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FIG. 1: (a) DOS, (b) MT, (c) AL , and (d) twisted AL dia- 
grams, plus (e) series of modified AL diagrams, that are used 
to calculate the compressibility and spin susceptibility. The 
four vector Q is taken to vanish in the static limit. 



interaction active above T c among correlated fermion 
pairs, in order to prevent the compressibility from di- 
verging when approaching T c similarly to what occurs 
for point-like bosons. To this end, it is necessary to 
improve on the standard t-matrix approach, which has 
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been successfully used in a variety of contexts but would 
now lead to a diverging compressibility when approaching 
T c from above. This is achieved by including the resid- 
ual interaction via the diagrammatic approach of Ref . Q 
which is equivalent to the Popov approach for (compos- 
ite) bosons in the BEC limit. The good agreement we 
will obtain by this approach with the experimental data 
on the compressibility (see below), over the whole tem- 
perature range where they are available, indicates that 
the residual interaction among correlated fermion pairs 
above T c Q represents a key ingredient for the thermo- 
dynamic stability of the system. 

For the spin susceptibility, only the DOS diagram of 
Figflja) and the MT diagram of FigQJb) are relevant 
above (as well as below) T c On physical grounds, for 
a Fermi gas with attractive interaction one expects the 
spin susceptibility to be suppressed when T c is reached 
from above and to vanish eventually deep in the super- 
fluid phase when T « T Cl by the argument that partners 
in Cooper pairs get locked in a spin singlet. This feature, 
which can be obtained already at the mean-field level in 
a purely BCS approach should be even more pro- 
nounced by the occurrence of pairing fluctuations above 
(as well as below) T c Q, which are associated with the 
occurrence of a pseudo-gap. This behavior results also 
from our calculation for the attractive gas when pairing 
fluctuations are included, but is not consistent with the 
experimental data for the spin susceptibility of Ref. [l| . 

In Ref.[l[ the lack of suppression of the spin suscepti- 
bility close to T c was seen as challenging the existence of 
a pseudo-gap in the unitary Fermi gas [7[. However, mea- 
surements that probe directly the single-particle excita- 
tio ns [1 , 0] as well as a number of theoretical calculations 
lave supported the existence of a pseudo-gap in 
the unitary Fermi gas. We shall show that the above 
apparent contradiction can be resolved by assuming that 
the measurement of the spin susceptibility of Rcf.[l| ac- 
tually explores a non-equilibrium state associated with 
a quasi- repulsive Fermi gas, and therefore cannot be di- 
rectly compared with equilibrium calculations for an at- 
tractive Fermi gas, as remarked already in Ref. 14 1. 

Although the experimental data for the compressibility 
and spin susceptibility were reported in Ref.[l| only at 
unitarity, we shall extend our calculations to both sides of 
the Fano-Fcshbach resonance for couplings (kpap)^ 1 > 
on the BEC side and (kpap)^ 1 < on the BCS side of 
the crossover. Here, kp is the Fermi wave vector related 
to the density by n = k F /(3ir 3 ) and ap the scattering 
length. For the attractive gas, this extension to both 
sides of the crossover is required to compare with the spin 
susceptibility data reported in Ref.[15[ for the trapped 
gas. For the quasi-repulsive gas, one needs to extend the 
calculation up to (kpap)^ 1 = 10 to recover the results of 
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a "dilute" repulsive gas [16| with good accuracy. 

Figure [5] shows the temperature dependence of the 
compressibility at unitarity. The experimental data from 
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FIG. 2: Compressibility at unitarity vs T/Tf, normalized 
by the value for an ideal Fermi gas at T = 0. The data 
from Ref. [l[ (circles) are compared with alternative theoreti- 
cal calculations. Full line: most complete calculation, includ- 
ing diagrams a+b+c+d+e of Fig. 1; dashed line: calculation 
including diagrams a+b+c+d; dashed-dotted line: virial ex- 
pansion; dotted line: non-interacting Fermi gas. The inset 
compares the results of our most complete calculation when 
(^f^f) -1 = +1.0 (full line) and (fcpfflF) 1 = —1.0 (dashed 
line) with those of a non-interacting Fermi gas (dotted line). 



Fig.4(a) of Ref. Q (circles) are compared with the theo- 
retical results obtained for an attractive Fermi gas from 
the static limit of the density correlation function Xn in 
the normal phase above T c (in this way, only a couple 
of experimental data at the lowest T are missed). The 
calculations neglect (dashed line) or include (full line) 
the residual interaction among pre-formed Cooper pairs 
above T c , and are based, respectively, on the DOS plus 
MT and AL diagrams and on the DOS plus MT and 
the whole series of AL diagrams of FigfTJ The results 
obtained from the high-temperature (virial) expansion 
of Ref.[I3| (dashed-dotted line) and those for a non- 
interacting Fermi gas (dotted line) are also shown for 
comparison. Note how the inclusion of the residual inter- 
action among pre-formed Cooper pairs is essential to get 
meaningful results for the compressibility, which would 
otherwise diverge when T — > T c + within the standard 
t-matrix approximation fl8j . 

Figure [3] shows the corresponding temperature depen- 
dence of the spin susceptibility at unitarity. The exper- 
imental data from Fig. 4(a) of Ref.[l[ (circles) are com- 
pared with the theoretical results obtained for an attrac- 
tive Fermi gas from the static limit of the spin corre- 
lation function Xs- The results shown by the full line 
are obtained above T c by summing the contributions of 
the DOS [FigOTa)] and MT [FigOJb)] diagrams, and ex- 
tended to the superfluid phase by adding the fluctuation 
contributions associated with these diagrams below T c to 
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FIG. 3: Spin susceptibility at unitarity vs T/Tf, normalized 
by the value for an ideal Fermi gas at T = 0. The data from 
Ref. PU (circles) are compared with alternative theoretical cal- 
culations for an attractive Fermi gas. Full line: calculation 
including diagrams a+b of Fig. 1 above T c and their exten- 
sion on top of the BCS contribution below T c ; dashed-dotted 
line: virial expansion; dotted line: non-interacting Fermi gas. 
The inset compares the data of Ref. [l5|] (circles) for a trapped 
Fermi gas at equilibrium with our calculations (diamonds). 



the BCS result @,0] (see also 0]). The results obtained 
from the high-temperature (virial) expansion of Ref. 17 1 
(dashed-dotted line) and those of a non-interacting Fermi 
gas (dotted line) are again shown for comparison [lij . 

Our results reproduce well the virial expansion for an 
attractive Fermi gas up to T « 5Tp, and are consistently 
below those for a non-interacting Fermi gas. This indi- 
cates a tendency toward pair formation in the normal 
state, that at lower temperature leads to a pronounced 
drop in the spin susceptibility which signals the presence 
of a "spin gap" well before the onset of the superfluid 
phase [6|. The spin susceptibility vanishes eventually for 
T — > 0, reflecting Cooper pairing in spin singlets 

Our results for the attractive Fermi gas, however, show 
marked deviations from the experimental data of Ref. [l[ , 
which lie instead above those for a non-interacting Fermi 
gas also at high temperature and do not perceive the 
expected suppression due to the singlet correlation in 
Cooper pairs at low temperature. This may indicate 
that the specific dynamical conditions through which the 
spin susceptibility data were determined in Ref.[l[ have 
de facto excluded the formation of the two-body bound 
state, resulting in the formation of a "quasi-repulsive" 
Fermi gas with an effective repulsive interaction [l4| . 

Before passing to describe such a quasi-repulsive Fermi 
gas to connect with the spin susceptibility data of Ref. [l[ , 
it may be relevant to compare our theoretical calculations 
for the the spin susceptibility of an attractive Fermi gas 



with an alternative set of data from Ref. [l5j , which was 
taken at thermodynamic equilibrium for a trapped gas 
across the BCS-BEC crossover. These data are reported 
in the inset of Figj3] (circles) and are in good agreement 
with our theoretical calculations for the trapped system 
(diamonds) [the two couplings on the left (right) corre- 
spond to a temperature 0.l3Tp (0.191V)]. 

A route to the description of a repulsive Fermi gas 
star ting from an attractive one was provided recently in 
Ref. [20( , by the exclusion of the bound-state contribution 
from the density equation on the BEC side of unitar- 
ity within a Nozieres-Schmitt-Rink (NSR) approach. It 
turns out, however, that this approach results in a wide 
forbidden region of the temperature-coupling phase dia- 
gram, in such a way that it cannot be used to obtain the 
spin susceptibility close to unitarity at low temperature. 

Here, we propose an alternative approach to describe 
the quasi-repulsive Fermi gas, which focuses directly on 
the (particle-particle) ladder propagator Tq entering the 
response diagrams of FigJT] and eliminates the effects of 
the bound state on the BEC side of unitarity as follows. 
Let us consider the spectral representation 



+ OC 



dw Imr^(q, w) 

7T — UJ 



(1) 



for the ladder propagator of the attractive Fermi gas. 
Here, q is the center-of-mass wave vector, f2„ = lixvT [y 
integer) a bosonic Matsubara frequency, and r^(q, uS) — 
r (q, ift u — >• OJ + irj) with r\ = + . To exclude the con- 
tribution of the two-body bound state on the BEC side 
of the resonance, we need to eliminate from ImT^q, tu) 
the delta-like (polar) contribution at the given q. This is 
done by starting the ^-integration in Eq.([T]) at the con- 
tinuum threshold given by w c (q) = q 2 /(4m) — 2/i, where 
m and \x are the fermionic mass and chemical potential, 
respectively. 

An analogous restriction on the frequency integration 
in the expression of the density is only what was required 
within the approach of Ref. [201] to get the thermody- 
namics of the quasi-repulsive Fermi gas. Our use of the 
spectral representation ([1}, however, requires us to take 
also into account an additional frequency-independent 
two-body term, which needs to be subtracted from the 
many-body ladder propagator in order to reproduce the 
correct behavior of a weakly repulsive Fermi gas when 
(fcjrOF) -1 ^ 1. This frequency-independent term can be 
inferred from the work of Ref. [2l| , and yields eventually 
the following expression for the ladder propagator of a 
quasi-repulsive Fermi gas: 



p+oo J 

rHp.q,iO0 = -/ - 



dw ImT^(q,a;) 8ir/(maF) 



ifl u — LU 



(2) 

where 2p is the incoming relative wave vector. Note that 
Eqs.([T]) and © coincide at unitarity. The occurrence 
of an energy-independent term is familiar, for instance, 
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in the dispersion relation for the forward scattering am- 
plitude in scattering theory [22[. We arc going to use 
the expression ^ to calculate the thermodynamics and 
the response functions of this out-of-equilibrium system, 
provided the density equation admits solutions [18l |. 
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FIG. 4: Spin susceptibility at unitarity vs T/Tf, normalized 
by the value for an ideal Fermi gas at T = 0. The data 
from Ref.[lj] (circles) are compared with alternative theoreti- 
cal calculations for a quasi-repulsive Fermi gas. Squares: our 
extrapolated values for a quasi-repulsive Fermi gas; dashed 
line: virial expansion; dotted line: non-interacting Fermi gas. 
The inset shows the coupling dependence of the spin suscep- 
tibility for a given T(— 0.22>) (full line), together with an 
extrapolation procedure based on two different fitting func- 
tions. 



In particular, we have found that the spin suscepti- 
bility at a given temperature as a function of coupling 
has the typical shape shown by the full line in the in- 
set of Fig|4] This shape coincides with that of a truly 
repulsive dilute Fermi gas [l6| when (kpap) 1 ^> 1, 
and departs only slightly from it even when the coupling 
gets reduced down to (kpap)^ 1 ~ 2 (this lower coupling 
turns out to be almost independent of temperature [23[ ) . 
The spin susceptibility starts then to drop when the cou- 
pling is lowered further toward unitarity, approaching a 
value that correspond to an attractive Fermi gas at the 
given temperature (a value that can be reached only for 
temperatures lying outside the forbidden region in the 
temperature-coupling diagram). 

Our assumption at this point is that, by extrapolat- 
ing the shape of the spin susceptibility before it drops at 
{UfCLf) 1 ~ 2 in the way shown by the broken lines in the 
inset of FigS] (corresponding to two different fitting func- 
tions), we should end up by reaching an excited configu- 
ration as if an avoided level crossing were present, with 
a dynamics determined by Landau-Zener processes [3]. 
[Similar ideas were discussed also in Refs.[2^, (2(| while 
analyzing the competing instabilities towards Stoner fer- 



romagnctism and pairing.] Correspondingly, we assume 
that the dynamics of the colliding clouds in the experi- 
ment of Ref . [l| activates a number of Landau-Zener pro- 
cesses, such that the data there reported for the spin 
susceptibility at unitarity as a function of temperature 
can be directly compared with our extrapolated values 
obtained by the above procedure. 

Figure 0] compares our extrapolated values for the 
quasi-repulsive Fermi gas at unitarity (squares) with the 
data from Fig. 4(a) of Ref.[l[ (circles) over an extended 
temperature range. The error bars attached to our ex- 
trapolated values derive from the statistical uncertain- 
ties produced by the different fitting functions (like in 
the inset). On physical grounds, these error bars may 
be thought of as associated with the underlying presence 
of a large number of Landau-Zener processes mentioned 
in Ref.[26[. Reported in the figure are also the results 
obtained by the high-temperature expansion of Ref. [13] 
where the bound-state contribution has been subtracted 
off (dashed line) and those of the non-interacting gas 
(dotted line). The good comparison between the ex- 
perimental data and our theoretical results supports our 
treatment of the quasi-repulsive Fermi g well as the 
assumptions about the underlying dynamical processes 
that result in the spin susceptibility data of Ref.[l[. It 
is further relevant to mention that by our approach to a 
quasi-repulsive Fermi gas there is no evidence for a Stoner 
instability toward a ferromagnetic phase, in accordance 
with a recent experimental finding [28|. 

In conclusion, we have reported theoretical calculations 
for the compressibility and spin susceptibility over a wide 
temperature and coupling range, and compared them 
with recent experimental data for Fermi gases. For the 
compressibility, the residual interaction among fcrmion 
pairs turned out to be an essential ingredient for compar- 
ing favorably with the experimental data for an attractive 
Fermi gas. The spin susceptibility was calculated both 
for an attractive and a suitably defined quasi-repulsive 
Fermi gas, and has been favorably compared with differ- 
ent sets of experimental data in the two cases. For the 
quasi-repulsive gas this comparison has relied on assum- 
ing a dynamics determined by Landau-Zener processes. 
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Supplemental material 
Diagrammatic formalism and linear response 

The number density 

P(r)=5>£(r)^a(r) (3) 



and spin density (h = 1) 

S*(r) = i5>£(r)aS^« 



(4) 



a,/3 



operators, where is a Pauli matrix, are expressed 
in terms of the fermionic field operator "0 Q (r) at space 
point r and with spin index a = (t, 4) ■ In terms of these 
operators, one writes the density-density 

X„(rr,rV) = -(T T [p(rr)p(r'r')]) (5) 

and spin-spin 

Xs (rr,rV) = — (T T [S z (rr)S*(rV)]> (6) 

correlation functions with imaginary time r. Here, T T is 
the time-ordering operator and 



T (7) 



is a modified Heisenberg picture with Hamiltonian H, 
number operator N, and chemical potential fi. Accord- 
ingly, in Eqs.([S]) and ([5]) the symbol (• ■ •) corresponds to 
a grand-canonical average The correlation functions 
([5]) and (j6]) representparticular cases of the two-particle 
correlation function [2|. 

From the above expressions, for a homogeneous system 
the (isothermal) compressibility Xn and spin susceptibil- 
ity \s follow as the static limits: 



dn 

Xn = T; 



dM 
~~b~h 



= lim x«(q, &u = 0) 

j, q-S-0 



lim x s (q,rJ„ = 0) 

q->0 



(8) 



(9) 



where n is the density, M the magnetization, h a uniform 
magnetic field, and 



nl/(fc B T) 



Xn/i (q,fi„) = / d{j - T >) e ^-r> (10) 

Jo 

x /"d(r-r')e-^( r - r ') Xn/s ( rr , r V). 

Here, q is a wave vector, Sl„ = 2ttvT (v integer) a bosonic 
Matsubara frequency, and ks the Boltzmann constant. 
The values of the static limits and © can conve- 
niently be normalized in terms of the corresponding non- 
interacting values 2iVo and 27Vo^, in the order, where 



A^o = mkp I (2ir) 2 is the density of states per spin com- 
ponent and fiB the Bohr magneton. 

Standard diagrammatic methods can be used to calcu- 
late the correlation functions x„(q, and Xs(q, ^i/) in 
Fourier space [l[ , both in the normal phase above T c and 
in the superfluid phase below T c . 

Choice of diagrams for 
compressibility and spin susceptibility 

Above T c , a pairing- fluctuation approach that ex- 
tends the Galitskii theory 0] throughout the BCS- 
BEC crossover Q identifies the relevant fermionic single- 
particle self-energy El with the diagram depicted in the 
upper panel of FigO where the particle-particle (ladder) 
propagator To is depicted in the lower panel of the same 
figure. It sums up all the elementary scattering processes 
between two fermions propagating in the medium with 
opposite spins owing to the contact nature of the inter- 
particle interaction. 



,(q) 



q-k q-k' 



-»- 






> 1 ) 




T (q) 















q-k q-k' q-k q-k' 

) i > t > 



k k' 



FIG. 5: Single-particle fermionic self-energy in the normal 
phase (upper panel) expressed in terms of the ladder propa- 
gator To between two fermions of opposite spins (lower panel). 
Full and dashed lines represent the fermionic propagator and 
interaction potential. 

The two-particle response that bears on this self-energy 
contains the effective two-particle interaction of the kinds 
depicted in FigJH] (that corresponds to Fig. 3 of Ref.[H[). 
To the lowest order, these terms produce the Aslamazov- 
Larkin (AL) diagram of Fig. 1(c) plus its twisted compan- 
ion of Fig. 1(d) of the main text, and the Maki-Thompson 
(MT) diagram of Fig. 1(b) of the main text. The two AL 
diagrams give equal contribution to the compressibility, 
but cancel each other for the spin susceptibility owing to 
the spin structure. On the other hand, the DOS diagram 
of Fig. 1(a) of the main text, where the above effective 
two-particle interaction does not appear, and the MT di- 
agram of Fig. 1(b) of the main text contribute to both 
quantities. This justifies the choice of diagrams made in 
the main text. 

Repeated structures based on the effective two-particle 
interactions of FigJS]arc also possible, and are specifically 
required by general conservation requirements [6[ . Along 
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cj > .... > a 



(a) 




(b) 



o < ■■■■ < o 



C > i i i * 
\ i i i o 



(c) 



FIG. 6: The three types of effective two-particle interaction 
obtained from the pairing self-energy of Fig[5j which give rise 
to the AL, twisted AL, and MT diagrams of Fig.l of the main 
text, in the order. 



these lines, for the calculation of the compressibility we 
have included the series of modified AL diagrams whose 
lowest-order contribution is depicted in Fig. 1(e) of the 
main text. In this case, these repeated processes are im- 
portant on physical grounds because they generalize to 
the BCS-BEC crossover analogous processes occurring 
for point-like bosons in the normal phase at the Hartree- 
Fock level, where they introduce the effects of the mutual 
repulsion between bosons and thus prevent the compress- 
ibility from diverging when approaching T c from above. 

In practice, in the static limit it is possible to in- 
clude this whole series of physical processes by exploiting 
a Ward identity that connects single- and two-particle 
fermionic Green's functions 0], whereby summing the 
whole series of modified AL diagrams is equivalent to 
calculating numerically dn/dfi in the following way. The 
density is obtained from the expression 

n = 2k B Tj2^ nn /(^3 G ( k ' w ») ( U ) 

in terms of the fermionic single-particle Green's function 
G(k, u„) = [ioJ n — k 2 /(2m) + ^ — Ep(k, w„)] _1 where k is 
a wave vector, u> n = irkBT(2n + l) (n integer) a fermionic 
Matsubara frequency, and 77 a positive infinitesimal. The 
fermionic self-energy Ep is formally of the type of El 
of FigEJ but with a dressed ladder propagator T that 
replaces To in order to include interaction processes be- 
tween composite bosons as described by the generalized 
Popov theory of Ref . Q . 
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FIG. 7: The diagrams of Fig.l of the main text are reproduced 
here with two thick lines each, which represent fermionic prop- 
agators dressed self-energy insertions. This dressing of the 
fermionic lines results when calculating the compressibility 
as \ n = dn/d/j, and the spin susceptibility as Xs = dM/dh 
through Ward identities. 



A comment is in order at this point about the degree 
of self-consistency that results when the diagrammatic 
structure for \n is generated in the above way by tak- 
ing the derivative of the expression (|11[) for n with re- 
spect to \x. One formally obtains the same diagrams of 
Fig.l of the main text that contribute to Xn (namely, di- 
agrams (a), (c), (d) and the whole series (e)), but now 
with two fermionic propagators (identified by the thick 
lines in FigJT]) which are dressed by the self-energy Ep. 

Below T c , the calculation of the spin-spin correlation 
function can initially be done at the level of the BCS 
(mean-field) approximation, whereby the bare bubble of 
Fig. 1(a) of the main text is replaced by the sum of two 
bubbles calculated, respectively, with two normal (G11) 
and two anomalous (G12) single-particle Green's func- 
tions In particular, in the static limit one obtains for 
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the spin susceptibility at the BCS level |1( 



(BCS) 



(T) 



1 

iVo 



dk df F (E k ) 



(2tt) 3 dE k 



(12) 



where S k = [(k 2 /(2m)-//) 2 + |A| 2 ] 1/2 i s the BCS disper- 
sion with gap A and Jf{E) = ( e E/{kBT) + the Fermi 
function. Note that this quantity vanishes in the zero- 
temperature limit, reflecting the singlet structure of the 
Cooper pairs. The mechanism for this to occur is a can- 
cellation of the contributions of the normal and anoma- 
lous BCS bubbles. 

Pairing fluctuations beyond mean field can be included 
below T c following the approach of Ref. ll|. In partic- 
ular, in the above two BCS (bubble) diagrams for the 
spin-spin correlation function the normal single-particle 
Green's functions Gn are affected by pairing fluctuations 
while the anomalous ones G12 remain at the BCS level. 
In addition, the MT diagram of Fig. 1(b) of the main text 
is introduced where now all single-particle lines are G\\. 
In this way one recovers the vanishing of the spin suscep- 
tibility at zero temperature for any coupling (kpap)^ 1 . 

Above T c , on the other hand, for the spin-spin cor- 
relation function the DOS diagram replaces the normal 
BCS bubble while the MT diagram plays the role of the 
anomalous BCS bubble. In this case, a complete can- 
cellation between the MT diagram and the fluctuation 
contributions to the DOS diagram occurs in the strong- 
coupling (BEC) limit for temperatures well below the 
pair-breaking temperature of the composite bosons. This 
is expected on physical grounds, since a non- vanishing 
contribution to the spin response for spin-less compos- 
ite bosons should result only when the temperature is 
comparable with their binding energy and the composite 
bosons break apart [Bj]. 

A comment on the degree of self-consistency for the 
diagrammatic structure is relevant also for \s- At the 



BCS level, no difference is introduced when calculating 
Xs via dMj dh with respect to the diagrammatic calcula- 
tion resulting in the expression (fT2|) . When fluctuations 
are introduced, while both fcrmionic propagators of the 
DOS diagram of Fig. 1(a) of the main text are affected by 
pairing fluctuations, only the pair of fermionic propaga- 
tors on the left side of the MT diagram of Fig. 1(b) of the 
main text are dressed by pairing fluctuations through the 
self-energy El depicted in FigEJ both above and below 
T c . This corresponds to the occurrence of thick lines in 
panels (a) and (b) of FigJT] We have, however, verified 
numerically that the dressing [TJb) of the MT diagram 
does not affect our main physical results, namely, that 
Xs is strongly suppressed for T <^ T c and vanishes at 
T = 0, while it slowly decreases for increasing tempera- 
ture when T > Tp where it coincides with the results of 
the high-temperature (virial) expansion. 
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